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Abstract. A one dimensional XX spin chain of finite length coupled to reservoirs 
at both ends is solved exactly in terms of a matrix product state ansatz. An explicit 
reprsentation of matrices of fixed dimension 4 independent of the chain length is found. 
Expectations of all observables are evaluated, showing that all connected correlations, 
apart from the nearest neighbor z — z, are zero. 
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1. Introduction 

Simplification leading to understanding of essential features of physical systems is one 
of the leading principles in theoretical physics. Finding simple solutions to seemingly 
complicated models is one of the ways to approach this goal. In quantum physics 
the complexity of a system grows exponentially with the number of particles. Even if 
an exact solution is possible representing it in a compact way is nontrivial for many- 
body systems. One approach to represent a given state is to express its expansion 
coefficients in a suitable basis in terms of product of matrices, the so-called matrix 
product states [21 H], used before in a wholly different context of 2d classical lattice 
model [3]. If quantum state is only weakly correlated the resulting matrices can be 
small. In condensed matter such an ansatz has been used to exactly describe ground 
states of many low dimensional spin systems, among the first one for instance the 
AKLT chain [4] or ladder system [5]. In addition to ground states, description in 
terms of products of matrices is also successfully used in algorithms for simulation of 
quantum systems [6]. Crucial for the efficiency of such simulation is the necessary matrix 
dimension D. Unfortunately, generic coherent quantum evolution will cause D to grow 
exponentially with the simulation time [7], rendering simulation inefficient. However, 
in some special cases of Heisenberg evolution of certain operators in integrable systems, 
like transverse Ising chain, the dimension D is small and does not grow with time [HE]. 
Efficient matrix product description with time-independent D is possible also for certain 
open versions of the same integrable models [9]. 

Besides quantum systems, matrix product states are also widely used in classical 
stochastic models. There a matrix product formulation is used to study nonequilibrium 
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stochastic lattice models, in particular their nonequilibrium stationary state (NESS). 
Since its first use (10] for the exact solution of a one- dimensional asymmetric exclusion 
process [Tl] a matrix product formulation has been utilized in a plethora of models, for 
a review see [12]. Expectations in the NESS can sometimes be calculated with the help 
of algebra only, sometimes though an explicit representation of matrices is required. It 
is not known in general when a finite dimensional matrix representation is possible. 

In the present work we provide an explicit finite dimensional matrix product 
solution for a stationary state of quantum system in nonequilibrium situation. This 
extends the applicability of matrix product states in quantum domain from ground states 
and time evolutions to the NESS. The model we consider is a one-dimensional XX spin 
chain coupled to reservoirs at chain ends. Using the Jordan- Wigner transformation one 
can reformulate the system in terms of spinless fermions. In the fermionic language the 
system is composed of free fermions. Hamiltonian part is therefore trivially integrable. 
In fact, for our choice of bath Lindblad operators the whole superoperator £ is quadratic 
in fermionic operators and can therefore be diagonalized in the operator space [Mj. 
An explicit solution for the first two orders in the driving can be also obtained as a 
particular case of an analytic solution for the XX chain with dephasing [14] . Slightly 
different model resulting in the same solution, in which XX chain repeatedly interacts 
with independent spins of the bath, has been explicitly solved in [15]. For studies of 
nonequilibrium states in a doubly infinite chain see [T71 H~6l ITS] . Our primary goal here 
is therefore not to solve the system but instead to provide a compact form for the exact 
NESS in terms of a matrix product ansatz with the matrices of low dimension D = 4 
that is independent of the chain length. This in turn enables one to explicitly evaluate 
the expectation value in the NESS of an arbitrary observable. 

Hamiltonian of the XX spin chain is given by 

n— 1 
3=1 

with standard Pauli matrices; lower indices running from j = l,...,n denote a site 
position. Dynamics of the spin chain coupled to the environment will be described in 
an effective way using the Lindblad master equation [19] 

f f p = i[p, H] + Cr\p) + Cr\p) = C(p). (2) 

To induce a nonequilibrium situation we couple the chain at the first and the last site to 
a "bath" , modelled here by linear operator £ bath which is expressed in terms of Lindblad 
operators L\ 2 acting on the first site, and L^ 2 acting on the last site, 

£\t (p) = E ([Ll' R P, L L ™] + [^' R , pL L k ^}) . (3) 

k 

We take the simplest Lindblad operators of the form 

L\ = ^r L (l - | + p) at, L L 2=] Jr L (l + ^-p) trf, (4) 
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on the 1st site, while on the nth site we have 



Lf = ^r R (l + | + p) a+, Lf = ^r R (l - I - p) a~, (5) 

af = (crj d= i crj) /2. Such bath operators induce an imbalance in magnetization, causing 
a flow of magnetization from one to the other end. For bath operators only, the 
stationary state, for which C^ th (p) = holds, is diagonal in the eigenbasis of af, with 
the expectation value (af) = tr (pal) = p—%, and correspondingly for £ R ath . Therefore, 
parameter p, plays the role of average magnetization, p its difference between the right 
and left ends, while IY r are coupling strengths. Matrix representation of superoperators 
£ L a R h can be found in the Appendix A We are interested in a stationary solution of 
the whole master equation, C(p) = 0, which we call a nonequilibrium stationary state, 
NESS for short, and denote simply by p. 

2. Matrix product ansatz 

We shall write the NESS p with the matrix product operator (MPO) ansatz, 

p = Yn £ (i|4 ai) 4 Q2) ---4 an) |i)^r^ 2 (6) 

z ai,a2,...,a n 

Indices run over labels of Pauli matrices forming an operator basis, G {x, y, z, 1}, 
with the convention aj = tj, matrices A[ a ^ are of dimension D x D, while |1) is a 
-D- dimensional unit vector. We arbitrarily choose its components to be 5j^. 

The goal is to write the NESS in terms of as small matrices as possible. Our method 
of solution shall be the following. In Ref . [14] , an exact solution for an XX model with 
dephasing has been provided of which our current XX chain is a special limit (the limit 
of zero dephasing). It has been observed that in a model without dephasing there are 
no long-range correlations, that is, all connected correlations are zero. This leads us to 
think that one could perhaps construct an exact matrix product operator solution just by 
observing nontrivial one and two-point observables; in our case these are magnetization, 
current and z — z correlations. Therefore, based on the solution from [H], we are first 
going to construct an MPO that accounts only for few-point observables in the NESS. 
Then we are going to show that such MPO in fact provides an exact soluton, that is, all 
observables and not just few-point are reproduced correctly. To show this we are going 
to use the algebraic properties of the matrices found. 

In [H] it has been found that (for r L = T R = 1, p = 0) the solution is (upto 
normalization) 

7~l — 1 2 Tl — 1 

4 4 3=1 iD J=l 

In the above expression we write only linear terms in p and z — z term of order p? . All 
other are inessential for the following discussion^]- How can we write such state in terms 

% Two quadratic terms not written in eq.JT]) are ^(—^i + a n)Ylk=i( a k a k+i ~~ <J k a k+i) + 
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of products of matrices? As an easy overture, we start with a simpler operator p obtained 
by keeping in the NESS © only terms with ct|, i.e., dropping the current term in eq.©. 

( \ 2 

One can easily convince oneself that the matrices Af = %|1)(1| + |2}(1| - ^|1)(2| and 



.4 



(i) 



|1)(1|, where a\ = — 



i • 



% «2, 



,n-l 



Obuik = 0, while trivially A) 



(x,y) 



0. 



gives the wanted state. Therefore, for this simple operator MPO ansatz with dimension 
D = 2 would suffice. To describe also the current though, we have to allow for at least 
two additional basis states in matrices Aj. We found that matrices of size D = 4 are 
sufficient. To correctly describe all terms of order /i and fi 2 (also those of order that 
are not explicitly written out in eq.©) the following set of matrices gives the correct 
operator, 



A 



-t 2 







/ CLi 

1 



V o 

-p, 

( o 



1 

V o 

where we use the notation A( x ' y ) 



o \ 




o J 



p p p -p 










t \ 




o J 



A) 



A (y) 



R 



(J) 



/ 1 



V o 



o \ 



t 

o J 



v — R, R, R, — R, — R, 



( 



V i 



o \ 




o J 



with period 4, A^ 4 = A\ KJ , and can be concisely written as A 3 
while A (y) 

(JS]) and appropriately chosen three parameters t, a\, a n and abuik = a>2,...,n-i is an exact 
NESS solution for the Lindblad equation (j2]) with arbitrary Tl,r, /i, fi- To show this we 
are going to use an algebraic approach similar to the one used in solutions of classical 
stochastic processes. 

/TT\ 

Let us denote by ^ the superoperator from the commutator part of master 
equation corresponding to the nearest-neighbor 2-spin term in the Hamiltonian, 
o~faf +l + o~Jo~J + i in our case. Writing four MPO matrices in a vector of matrices as 
Ai = (A- , A[ y \ A\ z \ A^), we can form all 16 different products of two matrices at 
consecutive sites through Aj <g) Aj+i (a 4 x 4 matrix, each element being a product of 
two D x D matrices). Depending on the action of Cf^ +l on products of operators, i.e., 
on £P]_ x (Aj Cg) A i+1 ), see also the Appendix A , solving for NESS is relatively simple 
in two cases 



( A ( x ' y ) A ( x ' y ) 



Matrices A\ x) and A\ y) 



(8) 



are periodic 



cos f j - sin \j) P, 



— (cos ^j + sin |j) R. We shall show that the MPO with the above matrices 



i) if the NESS is separable or has only two-particle entanglement, like 
for instance in valence bond states, then we can have Cf^ +l (Ai <g> A i+1 ) = 0; (ii) other 
relatively simple situation is when £P]_ x (Aj <g) A i+ i) results in a divergence-like term, 
that is 



Ai <8) M i+ i 



with some matrices M,- 



(Mf\M^,Mr,MD. 



r(») 



Mi 
r(i)^ 



Ai+u (9) 
Note that this ansatz is a trivial 



inhomogeneous extension of a standard procedure used in classical nonequilibrium 



A matrix product solution for a nonequilibrium steady state of an XX chain 



5 



systems [2U]. The reason to allow for spatially dependent matrices is that we want 
to find the MPO solution with the smallest D. We find that representation with D = 4 
is possible irrespective of the chain length n [§]. Any inhomogeneous solution can be 
written as a site-independent one but with larger matrices. One possibility is to make 
block site-independent matrix A of size D ■ n out of site-dependent matrices Ai as 



/ 



.4 








An 












A 2 









A n -! 





(10) 



If equation ([9]) holds terms from consecutive C^f +l will pairwise cancel, except at the 
boundaries. To ensure stationarity we have to enforce an additional condition at the 
boundaries. For our case of bath acting on a single boundary spin we get two equations, 



(1| 



Cl ath (A 1 )-M 1 =0, 



C h R ath (A n ) + M n |1) = 0. 



(11) 



If one manages to find such Ai and the associated Mj that eqs.([9J) and (fill) are satisfied 
one has found the NESS solution of the master equation (T5]). The condition (Q can 
in our case of the Pauli basis of the operator space be written as a set of 4 2 matrix 
equations of size D x D. As we have only 2 • 4 matrices Ai and Mj of size D x D, with 
8D 2 unknown parameters, it is not guaranteed that the solution exists. 

We are now going to show that for the XX model it actually does exist. We are 
going to find an explicit representation of matrices Mj, showing that they, together with 



Ai ((HI), satisfy equations (jHJ) and (ITTi) . Similarly as for A\ x ' 7 \ matrices M^ ,J> also have 
the periodicity 4. Their explicit (D = 4)-dimensional representation in the bulk, that is 



r( x >y) 



for the sites i = 2, 



n — 1, is 



(i) 



M 



M (x) 



S 



-2a 



bulk 



(S,S,- 

( o 




V o 



/ o 




V o 












o \ 



1 

o J 



-1 








M, 



M (y) 



T 



( 2t \ 
2t 
1 

V o o 10/ 

(T,-T,-T,T,T,.. 
( \ 



/ 






V o 



-1 









o / 



(12) 



One can convince oneself by direct calculation that the above Mj together with Ai (JH]) 
satisfy condition ([9]), written out in full the Appendix A eq. (1A.4I) . To also fulfill the 
two boundary conditions fTTTl) . M\ and M n must have additional matrix elements and 
parameters t, a±, a n and abuik must take specific values. We have M[ = S + 2(abu 



ai)|l)(4|, M[ y) =T + 2(a bulk -a 1 )|l)(3| and M x (z) = M 2 W + 4t(a bulk - a 1 )|l)(2|. On the 
§ For an explicit small dimensional MPO construction of some simple operators see [21] . 
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right end we have = M 2 (z) -f (a bulk -a n ) |2) (1|, M« = sgnx n (S+f (a bulk -a n )|3)(l|), 
AfW = sgny n (T + |(a bu i k — a„)|4)(l|), where sgnx n and sgny n are the signs in front of 
last S or T in eq.( !T2|) and depend on the site n. In addition, for the boundary terms to 
be zero the values of parameters must be 

r L r R 



t — ix- 

ai = fx 



i + r L r R )(r L + r R )' 
n (r L - r R ) + r L r R (r L + r R ) 
2 (i + r L r R )(r L + r 



R) 



0'2,...,n-l — Obulk — M — 7T 



n (r L -r R )(i-r L r 



2 (i + r L r R )(r L + r R ) 
_ fj, (iy- r R ) - r L r R (r L + r R ) 
"^2 (i + r L r R )(r L + r R ) ■ {U) 

Note that we have a\ — abuik = — £Tl, a n — abuik = ^r R and a n — a,\ = t(I\ + T R ). 
Parameters (ITS]) together with matrices (|S]) form an exact MPO solution of the NESS 
for the XX chain. 

If we add a homogeneous magnetic field in the z-direction to our Hamiltonian, 
that is the term of the form BEiLi '?! the model can again be solved exactly. In 



1 1 



fact, the single-site superoperator due to the magnetic field acts as C\ (erf) = 2B a 
Cf\al) = -2Baf, while ^(o? 1 ) = 0. Due to the symmetry of the NESS without 
the field (see also the explicit form of all nonzero terms given in [H]) and the minus sign 
in the action of C[ B \ we have T,i^T\p) = 0- This means that the solution presented 
OH]) is also the exact NESS in the presence of an arbitrary homogeneous field of strength 
B. 



3. Expectations of observables 

With an explicit representation of matrix product solution at hand we can evaluate 
expectations of various operators in the NESS. Let us first evaluate expectation of a 
series of o\ operators, like (otcr^ • • •). One or two er z s are simple to evaluate by direct 
calculation, we obtain (of) = and (ofcrj) = a^aj — t 2 5 i+ ij. Connected correlation 
function of n operators, (0^0^ ■ ■ ■ Oi n ) c , is obtained from an ordinary expectation 
value by subtracting product of all connected correlations where each involves less than 
n operators. For instance, (Oj) c = (Oj), (OiOj) c = {OiOf) — (Oi) c (Oj) c , (O^O^O^c = 

so on. Connected correlation function of two <r z s is therefore nonzero only on neighboring 
sites, (ofapc = — t 2 5i + \j. We are now going to show by induction that all higher order 
connected correlation functions involving more than two a z s are identically zero. Assume 
that the statement holds for products of upto n <r z s. Denoting by Z = a z cr z • • • <y\ = Za\ 
a product of n not-necessarily neighboring <r z , we would like to show that the connected 
correlation (Zo~f) c is zero for n > 2. Assuming connected correlations of more than two 
cr z s are zero, we can write 
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= (Za?)-(Z)(a?) c -(Z)(ala?) c , (14) 

where we denoted by Yf(Z) c a sum of all products of connected correlations, each 
involving < n operators <r z , for instance, ^'(of cr|) c = (cr 4 z cr|) c + {af) c {a^) c . We used the 
fact that to have a nonzero connected correlation o~f must be paired with at most one 
other o~j and that it must be its neighbor (non nearest-neighbor connected correlations 
are zero). Using our explicit MPO representation (jHJ) of matrices for and A^ we 
will now show that the right-hand side of (fT4|) is zero. First, observe that the matrix 
corresponding to the product of matrices occurring in the operator Z has an upper-left 
2x2 block equal to 

„ ( (z) -e(z) \ 



(15) 



Upper left element is (Z) by definition while the upper right follows through a simple 
multiplication of matrix corresponding to Z by A k . With the explicit form of Z we 
have 

(Za?) = (Z)(o-t) c -t 2 (Z)5 k+hh (16) 
where the Kronecker delta takes into account multiplication of the matrix for Z by A^h 



in the case of non-neighboring sites k and I. Plugging this into eq. (|14p we see that the 
(n + l)-point connected correlation of cr z is indeed zero. This completes the proof. 

Single point expectations (erf) and (erf) are zero. Among two-point expectations on 
neighboring sites, besides z — z, the only nonzero terms are two from the spin current, 
that is (o~fcrJ +1 ) = —(o~Jo~f +l ) = —t. Therefore, the expectation value of the spin current 
operator j k = 2(cr£of +1 —a y k a^ +l ) is (jk) = —At. If we have a product of non-overlapping 
current operators at different sites, (jijkjr 4 -) (with k > i + 1,1 > k + 1, . . .), the 
expectation value is simple. Observe that A^A!$ 1 = — £(|1)(1| + |3)(3|). Because of 
the form of A^\ if we have a product of only A^s and p terms A^A^ 1 at various 
sites i, their expectation value is simply equal to t v . This means that the connected 
correlation function of non-overlapping operators ji is zero, (jijkji • • -)o = 0, apart 
from (ji) c = —At. If two current operators overlap, i.e. jiji+i, with the hermitean 
(jiji+i + ji+iji)/2 = —A(afo-f +2 + °i°i+2)i ^ ne corresponding product of matrices is 
A!p A$ x a!$2 = — 1 1 > < 1 1 +t|3)(4|). From this form, similarly as before, we can 
see that for products of ji the term |3)(4| is not important, resulting again in all 
connected correlations being zero. Similar argument holds for products of more than 
two overlapping current operators. 

Finally, let us discuss expectations of products of a\ and jk, the only remaining 
nonzero terms in the NESS. For overlapping sites, jjcrf + afji = (the same for , 
and we have to consider only non-overlapping operators. Two-point expectation is 
{ji a j) = —Ataj, connected correlation being therefore 0. Again, taking into account 
that if we consider only products of Aj, A^ and A^A^, the term |3)(3| in A^A^ t 
is irrelevant because there are no |3)(1|, |4)(1| or their hermitean conjugates in any of 
the three matrices. The expectations are therefore trivial products of scalar quantities 
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in front of |1)(1| terms, resulting in all connected correlations being zero. If we 
have a product of neighboring currents, like for instance in (jiji+i + ji+iji)/2, similar 
argument holds. Incidentally, we also see that the expectation of the energy current 
jf = 2«_iO r M+i - °i-i° z i°f+x), such that i[afaf +1 + crjaj +1} H] = jf - jf +l , is zero. In 
the NESS with our choice of baths no energy current flows. 

Another way to write the NESS is to rewrite it in terms of an exponential function 
as p = exp (-H). Doing the calculation we observe that the operator H contains only 
few nonzero terms. These are o~\ at all n sites (with different prefactors on different 
sites), spin current jk = 2b k at all n — 1 sites (the same prefactor on all sites) as well as 
hf k+1 \ with k = 1, . . ., and bf ^ with k — 2, . . ., where hf = o-Jo-j +1 ■ ■ ■ o'j +k _ 2 aJ +k _ l + 

^J^i+i ' ' ' cr i+fc-2 cr J+fc-i an d bj ^ = o-jO-j +1 ■ ■ ■ o"| +fc _ 2 o"J +fc _ 1 — a^a^ +l ■ ■ ■ o-j +k _ 2 o~j +k _ 1 . 
This shows that the NESS in the open XX chain can not be exactly written as 
a quasi-equilibrium generalized Gibbs state p ~ exp (J2j K jQj) with locally varying 
fields Kj, and Qj conserved quantities of the corresponding Hamiltonian system [22J. 
Namely, conserved quantities in integrable systems can depend on small changes in 
the system, for instance on boundary conditions. For Id XX chain with periodic 
boundary conditions two infinite sequences of conserved quantities exist, one set are 
Qtk = £j hf k) and Qt k+1 = Ej bf k+1 \ another Q 2k = Ej bf k) and Q 2k+l = £j hf k+1) . 
If we change boundary conditions to open, only half of conserved quantities survive, 
that is, for open boundary conditions conserved quantities are only Q 2 k = Qtk + l^k 
and Q2k+i = Q^+i + <?2fc+i, where qj are additional boundary terms [23J. First few 
conserved quantities for open XX chain are Q2 = Y%=1 hf , Q3 = Y^fZi hf + o\ + cr^, 
Qi = E"=i hf - hf ] - h^, Q 5 = E"=i 4 hf + o\ - hf + al_ x - hf_ 2 . We can see that 
H can not be written as a sum of conserved quantities Q of an open chain neither that 
of a periodic one. A change in the boundary condition can therefore globally influence 
constants of motion of Hamiltonian system. This sensitivity also translates to open 
systems described by master equation: a change in boundary operators of the bath 
can have a large influence on the NESS. In short, for open integrable systems there is 
an ambiguity which "conserved" quantities Qj should one use in the generalized Gibbs 
ensemble, or, in other words, the Qj depend on the bath. Such non-universality, where 
the functional form of NESS does not depend only on the Hamiltonian but also on 
reservoirs, is probably generic situation for open versions of integrable systems. 

To summarize, we have shown that all connected correlations, apart from 
neighboring (alal +l ) C) are identically zero. This should not come as a surprise. In 
fact, working in fermionic language the system is quadratic, therefore all expectations 
of products of fermionic operators can be evaluated in terms of two-point expectations 
using Wick's theorem. 



4. Conclusion 



We have found an exact solution of an open XX chain in terms of the matrix product 
ansatz. An explicit (D = 4)-dimensional representation of matrices is found, enabling 
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to evaluate arbitrary expectations. All connected correlation functions, apart from 
the correlation function of magnetization at nearest-neighbor sites, are identically zero. 
The results presented extend the applicability of matrix product states to quantum 
nonequilibrium systems. Furthermore, the method of solution is an algebraic one 
borrowed from the field of classical stochastic processes where it has been used very 
successfully. It is hopped that this will lead to new exactly solvable nonequilibrium 
quantum systems. One such instance is an open XX chain with dephasing, a solvable 
diffusive model [H] where, based on our experience, a compact matrix product solution 
is also possible. Unfortunately though, in an even more interesting XXZ model the 
algebra seems to be more complicated. 
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Appendix A. 

Appendix A.l. Representation of bath super operators 

Using a basis of Pauli matrices, and tensor products thereof, a single-site superoperator 
for the bath is 

/ -2 \ 



C 



bath 
L 



(A.l) 



0-2 
0-4 -2(//-2/2) 
V J 

Superoperator £^ ath for the right bath is obtained by replacing Tl with Tr and \x by 
Basis elements are ordered as (<r x , <r y , a z , 1). Matrix representation gives us the 
operation of £ bath on the expansion coefficients of operators in the Pauli basis, e.g., 
writing p = J2 c aO~ a , we have Cp = p' = Z)c^o" Q , with c' a = J2p £<x,p c /3- Alternatively, 
the same matrices can be thought of as transforming matrices in the MPO ansatz, e.g., 

Appendix A. 2. Representation of C ii+1 

Superoperator transforms operators according to their commutator with the 

Hamiltonian. It is therefore fully specified by its operation on a 16 dimensional basis of 
2-site operators. We choose products of Pauli matrices as a basis. As an example, for 
the XX chain we have for instance, £-^]_ 1 (ljO"f +1 ) = — 2<xfof +1 . Acting on MPO ansatz, 
this for instance gives • • A^A^ ■ ■ •) haf +1 ) = (■ • • {-lA^A^) ■ ■ •) afa^ v 

In the transformed state, the term that comes in the product of matrices at the position 
of At Afh is therefore equal to (— 2A ( i 1 ' 1 A^ x ). Because we are interested in how Cf^+i 
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transforms MPO, it is handy to represent its action in terms of a 4 x 4 matrix whose 

as 

/ 



element [-C^li]^ gives us the term in front of the c"cr? +1 . For the XX chain we have 



/•(H) 



where B = A^A+i - Af> Af>^ and C 
transformation 





-B 

A i ^i+l 

Ay) 4(2) 

|(1) aV) 



B 


_ /l( x ) a( z ) 



-A w A [x) 



_a(z) 4W 





C 



-C 




(A.2) 



/ 



4W A (y) 



(■ ■ ■ ([A,i+i]tti,«t+i) ■ ") a i t<7 i+i 



A^A^i, explicitly giving the 

(A.3) 



Written explicitly, the algebra fl9]) of matrices Aj and Mj in the bulk that is induced 

'(H) 



by C\i+i is for the XX model 

2(4°^! 

-2(A (Z) 4S 

iW A (y) 



2(A 
-2(A. 



W 4 (y) 



i+l 

(y) 

i+l 



^i ^i+1, 
4(1) ^ 

^(y)^W n 



A- y) A 



i+i. 

(x) > 
i+l/ 









z/1 i ^i+i 

o 4(2)4 (y) 
z/1 i ^i+l 

-2A (1) A ( + \ 

_9 4 (z) 4 (x) 
z/1 i A i+i 

9 4 (y) 4(1) 
z/1 i ^i+i 
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Representation given in eq.® and (IT21 satisfies this algebra for the arbitrary values of 
four parameters ai, a n , Obuik and t. 
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